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Abstract 

Effective action is proposed for the problem of Coulomb blocking of tunneling. 
The approach is well suited to deal with the "strong coupling" situation near 
zero bias, where perturbation theory diverges. By a semiclassical treatment, 
we reduce the physics to that of electrodynamics in imaginary time, and 
express the anomaly through exact conductivity of the system a(u, q) and 
exact interaction. For the diffusive anomaly, we compare the result with the 
perturbation theory of Altshuler, Aronov, and Lee. For the metal-insulator 
transition we derive exact relation of the anomaly and critical exponent of 
conductivity. 
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Introduction. Coulomb correlations in a disordered metal suppress tunneling conductiv- 
ity at low bias, and lead to the "zero-bias anomaly" that has been known experimentaly 
since the early sixties [[I]]. A perturbation theory of this effect was developed by Altshuler, 
Aronov and Lee 0. The theory deals with the diffusive limit, and shows that the blocking 
increases at small bias which leads to a singularity in the tunneling conductivity. The theory 
has been thoroughly tested experimentally and found to be extremely accurate when the 
anomaly is a weak feature on the top of a large constant conductivity. 

However, in the past few years there has been interest in the systems with strong Coulomb 
effects, such as disordered metals and semiconductors near metal-insulator transition ||. It 
was found that the Coulomb anomaly sharply increases near the transition, and thus it 
provides a test of the role of Coulomb correlations in strongly disordered systems. Another 
important discovery is the observation by Ashoori et al. of the Coulomb blocking of tunnel- 
ing in a two-dimensional metal in magnetic field [3J. In this experiment it was found that 
at certain magnetic field the zero-bias anomaly abruptly increases and transforms to a "soft 
Coulomb gap" . It has been pointed out H that this transition is induced by disorder. More 
recently, the gap was studied in the systems with higher mobility [§,0, where current is al- 
most entirely blocked below certain threshold bias. These findings caused a lot of theoretical 
work, concerned with the behavior of current near zero bias [§], and with determining the 
gap width ||. 

It is characteristic for all works cited above that they start with some model of the 
system, and then proceed with a many-body calculation. However, it would be desirable to 
have an approach that does not rely on any model, which would allow a comparison with 
experiment in the situations where there is no accepted model for conductivity. The goal 
of this paper is to propose an effective action theory that treats the Coulomb anomaly as 
cooperative tunneling and gives tunneling rate in terms of the actual conductivity a(u),q) 
of the system. For the first time, it provides non-perturbative treatment of the anomaly, 
accurate in the strong coupling regime. For example, at the metal-insulator transition, 
conductiviy has scale invariant form ct(uj) ~ u a , and we are able to predict the form of the 
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tunneling I — V curve. This relation may be used to determine the critical exponent a. 
We discuss it below, and also apply the method to diffusive anomaly in a two dimensional 
metal, with and without magnetic field. Comparison with the perturbation theory results 
is given, and an agreement is found. 

The physics of tunneling involves motion of a large number of electrons: while only 
one electron is actually transferred across the barrier, many other electrons are moving 
coherently to accomodate the new electron, and this collective effect completely controls the 
tunneling rate. We will argue that this motion can be treated semiclassically as classical 
electrodynamics in imaginary time, find instanton solution, and derive an expression for 
the tunneling rate in terms of the instanton action. To justify the semiclassical picture, let 
us consider a situation when at small bias one electron crosses the barrier. Typically, the 
time of flight through the barrier is much shorter than the relaxation time in the electron 
liquid. Therefore, while one electron is traversing the barrier other electrons practically do 
not move. Thus instantly a large electrostatic potential is formed, both due to the tunneling 
electron itself, and due to the screening hole left behind. The jump in electrostatic energy 
by an amount much bigger than the bias eV means that right after the one electron transfer 
we find the system in a classically forbidden state "under" the Coulomb barrier. In order to 
accomplish tunneling, the charge yet has to spread over a large area, so that the potential 
of the charge fluctuation is reduced below eV. If the conductivity is finite, the spreading 
over large distance takes long time, and thus the action of the whole system under-barrier 
motion is much bigger than h. 

An attempt to derive the Coulomb anomaly semiclassically was made by Spivak |7j 
who considered a two dimensional metal, and used Maxwell's theory of charge spreading, 
r(t) ~ at, to calculate the Coulomb part of the action: 




From that, the tunneling is suppressed by exp (— ~<S(i*)J, where £* = %/eV. The estimate 
shows that the action grows at small bias, which is a clear sign that the semiclassical treat- 
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ment is meaningful even for a well conducting metal. However, in the diffusive limit the 
estimate ([I]) does not agree with the perturbation theory. We shall see that the reason is 
that the main part of the action is rather Ohmic than Coulomb, and that after writing the 
action properly the semiclassical method completely recovers the pertubation theory result. 
The action. Let us write the action in terms of the charge and current densities p(r, t) 
and j(r, t). The contribution to the action of the spreading charge is mainly coming from 
long times when the charge deviation from equilibrium is small. Therefore, we can expand 
the action in powers of p(r, t) and j(r, t), and keep only quadratic terms. The action should 
reproduce the classical electrodynamics equations: the Ohm's law and charge continuity. In 
principle, this requirement is sufficient to determine the form of the action. However, it is 
more convenient to argue in the following way. We are going to use the action to study the 
dynamics in imaginary time. Therefore, the action is precisely the one that appears in the 
quantum partition function. The latter action expanded up to quadratic terms in charge and 
current density must yield correct Nyquist spectrum of current fluctuations in equilibrium: 

Here 

g=i + DVp (3) 

is external current and D a p is the tensor of diffusion constants related to the conductivity 
tensor by the Einstein's formula: a = e 2 uD, where v = dn/dp is compressibility. Generally, 
both a and D are functions of the frequency and momentum. For simplicity, we assume that 
the temperature is zero and discuss only a two dimensional metal with spatially isotropic 
and homogeneous conductivity: a xx = a yy , a xy = —<J yx . 

The requirement that the action produces correct current fluctuations is essentially equiv- 
alent to the fluctuation-dissipation theorem. Thus, in imaginary time we get 

S = - J J d A xid 4 X2 



Si K Xl _ X2 \ 



; 2 + 



(4) 
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where x 1)2 = (ti,2> ^1,2)5 an d $12 = 5(ti — t 2 ). The kernel K rt is related to the current 
correlator, 

(KiU = ((sl, q s-^-,)) ( 5 ) 

given by (0), where a and -D are functions of the Matsubara frequency related with the real 
frequency functions by the usual analytic continuation. We take Coulomb interaction in the 
second term of the action (|J) as non-retarded because we are going to study systems with 
relatively low conductivity, and thus slow charge relaxation. 

Dynamics in imaginary time. To calculate the tunnelling rate, we use the instanton 
method and look for a path in imaginary time. Among the "bounce" paths symmetric in 
time, p(r, t) = p(r, —t), j(r, t) = — j(r, —t), we shall find the least action path which will give 
a semiclassical estimate of the tunneling rate exponent. 

Let us derive equation of motion from the variational principle. We note that the action 
(H) contains the charge and current densities as independent variables, as Eq. @ was derived 
by matching with the equilibrium fluctuations in the grand canonical ensemble where charge 
is not conserved. Therefore, we have to supply the action (|4|) with the charge continuity 
constraint: p + V ■ j = J{r, t), where 

J(r,t) = e5(r)(5(t + r)-5(t-T)) . (6) 

This form of the charge source J7"(r, t) describes electron entering the system at r = 0, 
t = —t, and exiting at t = r through the same point. 
By incorporating the charge continuity we get 

Stotai = S(p,}) + 0(r, t) (p + Vj - J(r, t)) , (7) 

where 0(r, t) is Lagrange multiplier. For the least action path, the variation of Stotai relative 
to infinitesimal change 50, 8p, and <5j vanishes. (Note that due to the charge continuity 8p + 
V • (5j = 0.) After eliminating we get the standard equations of classical electrodynamics: 
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(i) p + V-j = J{r,t) ; 

(ii) } + DVp = a(uj,q)E; (8) 

(iii) E(r, t) = -V r J dr'p(r', t)U(\r - r'|) . 

The equations describe the system trajectory in imaginary time, i.e., the spreading of charge 
under the Coulomb barrier due to selfinteraction. 

The next step is to solve Eq.([|) for p and j, and to compute the action (Q). For a spatially 
homogeneous system, by using Fourier transform, we get 

J(u) 

PKUJA) \u\ + Dq* + a xx q*U q ' 
j(w,g) = -iK~ 1 (u,q)q L U q p(u,q) , 

where U q is the Coulomb potential formfactor. Substituted in Eq.([|) this yields the action 

S (r)-W |JM|2 U < (9) 

01 1 2^ \u\ + Dqi \u\ + Dq* + a xx q*U q {) 

which depends on the accomodation time r through Fourier component of the charge source: 

Jijjj) = 2ie sin out. 

To obtain total action of the system we subtract from the action Sq(t) of spreading charge 
the term 2eVr that accounts for the work done by voltage source: <S(r) = Sq(t) — 2eVr. 
Thus the energy conservation at transferring one electron across the barrier is assured. Then 
one has to optimize «S(r) in r. Optimal r* satisfies the relation 

g%g» = 2eV (10) 

OT 

Having solved Eq . (|10D for r*, one obtains tunneling rate that coinsides with conductivity up 
to a constant factor: 



G(V) = G exp 



■UsqMV))- -2eVr*(V)) 
n 



(11) 



The optimal t* can be interpreted as the charge accomodation time. 

The accuracy of the term 2eVr is determined by the assumption that ^> r/, the time 
it takes one electron to traverse the barrier. This assumption is valid whenever there is an 
anomaly: if t* ~ t/, then <S(r*) rs 7i, and thus there is no tunneling suppression. 
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Taken together, the equations @, fllOD , and ( pTTj) define conductivity. With this general 
framework, one can study the anomaly in different systems. Let us emphasize that after hav- 
ing calculated t*(V) and <S(r*) it is essential to check the selfconsistency of the assumption 
that t* ^> Tf. To illustrate this point, let us consider a clean metal, where u{oj) = ine 2 /muj. 
In this case Eq.(|9]) gives Sq(t) ^ h at any r ^> Tf, and henceforth r* ~ r/. Naturally, this 
indicates absence of the anomaly in a clean metal, the result familiar from the Fermi liquid 
picture. 

Diffusive anomaly. For a two dimensional metal with elastic scattering time tq and non- 
screened Coulomb interaction we set U q = 2ir/\q\ and a xx = a, constant at |u;|, wf|<7| < 1/tq. 
Then Eq.(|9|) gives 

5M =8^ ln ©' n ( TT » ffV2)2 ) ' (12) 
(Here v is compressibility.) From Eq. (|T0|) . 

n = ^ Hhave/V) . (13) 

The theory is selfconsistent in the hydrodynamic limit, > r , i.e., at eV < e 2 /cxT . Then 
the least action is 

e 2 / e \ / er a(z/e 2 '' 2 N 



^'^'Hsv^J'H^f^J (14) 

It is interesting to compare this result with the identical double-log dependence derived by 
Altshuler, Aronov, and Lee in a different context p|. They calculated perturbatively the 
correction to the tunneling density of states 5v{e) with the assumption that it is small, \5v\ 
z/o, which is the case only for a weak disorder. It was found that Su(e) = — h^v^SiV = 
e/e), where S{V) is given by ([14]). The main difference is that our double-log has to be 
exponentiated to get the tunnelling density of states, while in the double-log itself appears 
as a correction to the density of states. In the range of the perturbation theory validity the 
two results agree. From that point of view, our calculation provides description of the 
diffusive anomaly at low bias, where the perturbation theory diverges. 



Screening by electrodes. In a real experiment the charge tunnels between two elec- 
trodes, and there are separate contributions to the action due to the relaxation of the electron 
and hole charges on both sides of the barrier. If the electrodes are close, the charges partially 
screen the field of each other, which makes their spreading correlated. In this case the least 
action is smaller than the sum of independent contributions of the electrodes, and thus the 
anomaly is weakened. For a two dimensional system, in terms of algebra, the effect will be 
that the log- divergence of the integral over q in Eq.([|) will be cut at q ~ a -1 , where a is 
the distance between the electrodes. As a result, the V— dependence of the second log in 
Eq.flHD saturates at eV c^Vo = Tia/a. 

This "excitonic" correlation effect can be treated straightforwardly by writing the action 
(f|) for each electrode separately, together with the term describing interaction across the 
barrier. First, let us consider two identical parallel planes at distance a, and assume e 2 va ^> 
1, which is the case in almost all experiments. Then the least action is 

Sq(t) = a In (^j at r > K/eV . (15) 

Here a = ttt- In 2ne 2 i , a. If the planes have different conductivities and densities of states, 
the least action still has the form (|T5|), but now 

1 r 1 , Anaa 1 , Anaa' 



a 



An 2 



In — — 1 In 



(16) 



where a = 0102/(01 + 02), and the subscripts 1,2 label the planes. 

We find that at low V < Vq the I — V curve is given by the power law I ~ V a+1 . 
As expected, the tunneling suppression in this case is weaker than for the non-screened 
interaction. 

Anomaly near metal-insulator transition. The effective action theory remains valid 
even for a highly disordered system where the Drude model of conductivity does not work. To 
illustrate this point, let us consider disordered metal near metal- insulator transition. The 
problem has been extensively studied, and properties of conductivity can be summarized 
as follows [|l],|l(],|ri|] . Static conductivity of the critical state vanishes, and as function of 
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frequency it obeys scaling law: cr(uS) ~ u a , where the exponent a is constrained by (d — 
2)/d < a < 1. The lower bound a = (d — 2)/d is reached for a non-interacting system, and 
interaction shifts a up. It is convenient to write a = (d — 2)/(d — £), < £ < 2. 

We substitute the scaling form in the action (|9|) and get S(r) ~ r' 3 , where (3 = ((d — 
2)/2(d — (). Then the time t*(V) ~ V -1 /^-^ which leads to the tunneling rate 

G~e w -(V/ V y , l=^fy (17) 

with universal 7. The constant V depends on bare conductivity, and thus is not universal. 
Let us emphasize again that the conductivity of a strongly disordered metal is a difficult 
many-body problem which we do not attempt to address here. The meaning of the result is 
that it shows how conductivity, if known, can be used to analyze the anomaly. 
2D electron gas in magnetic field. It is straightforward to incorporate magnetic field 
in the theory by substituting a xx = a xx (B) in the action (0). (Note that a xy does not 
enter.) As magnetic field increases, the conductivity drops, and at certain field it reaches 
the quantum limit a q = e 2 /h. In this field range the prefactor a in Eq.(^) becomes of 
the order of one, and the anomaly in the conductivity changes from weak to strong. The 
threshold conductivity, according to Eq. (|T2"D , is 

1 e 2 

a c = — -— ln27re 2 z/a (18) 

A transition like that was observed by Ashoori et al. |4j] in the tunneling current from a 3D 
metal into a 2D electron gas. In this experiment, the ohmic conductance was measured as 
function of temperature, which corresponds to our zero temperature non-linear current taken 
at V ~ ksT/e. The 2D gas was relatively clean with the zero field mobility corresponding 
to the elastic scattering time To ~ 4 • 10~ 12 s. The Fermi energy calculated from the 
electron density was Ep ~ 10 mV. By using the result (|T8"| ) together with the Drude- 
Lorentz model, a xx (B) = ne 2 /mTQU 2 at Tquj c ^> 1, one finds that the anomaly hardening 
transition corresponds to the cyclotron frequency uj* = (SnEp/hro) 1 ^ 2 ~ 8.0 mV . In terms 
of the field intensity this is approximately 4.6 Tesla which is quite close to the transition 
field reported in Ref ||. 
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It is interesting to note that in a weakly disordered metal with EfTq 3> % the threshold 
field is small: TiuJq <C Ep. This means that the transition occurs well below the field where 
the Quantum Hall state is formed. Therefore, our estimate of Uq based on the "bare" 
Drude-Lorentz conductivity is meaningful and legitimate. On the other hand, to find the 
current at very low V one would have to use the conductivity renormalized by localization 
and interaction effects. 

Finally, let us mention a relation with the work by Halperin, He, and Platzman || that 
deals with the anomaly in the v = 1/2 Quantum Hall state. In this work, the problem was 
treated by summing linked cluster terms of perturbation theory, with the density responce 
function borrowed from the Chern-Simons Fermi liquid theory [jl2|| . The anomaly was found 
to have the form: 

G(V) ~ exp -Vq/V , V = Air-^i , (19) 

where V <S Vq, and n is density. It is interesting to see how this result can be derived from 
the effective action. It has been shown |l2| that conductivity of the v = 1/2 state has strong 
spatial dispersion: tx fc = A\k\, A = e 2 /IQne^/nn. If this form is inserted in the action @, 
one gets S(r) = tt^2t/A, which leads to the tunneling rate (1191) . 

Conclusion. We argued that the theory of the Coulomb anomaly in the regime of strong 
suppression of tunneling is semiclassical. The underlying reason is that the transfer of one 
electron across the barrier is controlled by cooperative motion of many other electrons. We 
treat this motion as classical electrodynamics in imaginary time, write the action and find 
instanton trajectory. Relation with the perturbation theory is discussed and an agreement is 
found. The tunneling current is expressed in terms of the actual conductivity of the system, 
which is useful in the situations where there is no accepted model of conductivity. 
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